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THEORETICAL METHOD FOR DERIVING M EARTH-CENTERED DATUM 


FROM OPTICAL OBSERVATIONS OF THE EARTH'S HORIZON 

FROM AN EARTH SATELLITE* 

By Ruben L. Jones 
Langley Research Center 


SUMMARY 


The necessary equations for a theoretical method for obtaining directly an 
earth-centered world geodetic system from angular measurements involving the 
earth's visible horizon as seen from an earth satellite are derived. As a 
check on the method, a computer program was written whereby data could be sim- 
ulated and results obtained. 

The method is found to be mathematically feasible, and the results for a 
particular case are given to show the accuracy which can be expected under 
ideal conditions. Further, when large variances in the observations were 
assumed, the method was found to converge, although slowly. Finally, some of 
the conditions for the most effective application of the theory are determined. 

Althou^ the method is found to be mathematically feasible and to converge 
for large variances in the measurements, the results are preliminary and a 
detailed statistical analysis is needed. As for additional applications, the 
theory can be useful in a variety of theoretical studies, such as a study of 
the accuracy penalties resulting from assuming the earth to be spherical rather 
than oblate in navigation schemes and horizon-uncertainty studies. 


INTRODUCTION 




The technological advances since World War II have instigated a phenom- 
enal growth in interest for an accurate world geodetic system. This increased 
interest has resulted chiefly from the national objective of space exploration. 

In order to complete a deep space probe, it is desirable to know the pre- 
cise locations of the launch site and of the radar stations which will observe 
the flight of the space vehicle. The error in the location of a point on the 


*The material presented herein is based on a thesis submitted in partial 
fulfillment of the requirements for the degree of Master of Science in Geodetic 
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earth* s surface is not necessarily a function of the error in the survey alone. 
Because of the irregular shape of the earth* s surface^ the variation of the 
earth* s internal density vith location, and the uncertainty in the location of 
the earth* s center, the accuracy with which the launch site and radar stations 
can he located on a particular datum is limited. 

Geodesists have developed several sophisticated methods for dealing with 
each of these factors. Even so, to date, data are insufficient to determine 
precisely the earth's center and the ellipsoid which will best fit the earth 
as a whole. Recently, the artificial earth satellite has been employed for 
geodetic purposes. A major finding has been that the flattening f is approxi- 
mately 1/298.3 rather than the 1/297 previously supposed. Further, the last 
three digits of the magnitude of the semimajor axis a, previously thought to 
be 6 378 388 meters, are seriously questioned. 

The necessary equations for a theoretical method for obtaining directly an 
earth-centered world geodetic system will be derived herein. The results for 
a particular case will be given to show the accuracy which can be expected 
under ideal conditions. Further, it will be determined if the method will, in 
fact, converge when large variances in the observations are assixmed; finally, 
some of the conditions for the most effective application of the theory will 
be determined. 


SYMBOLS 


A,B,C coefficients of quadratic equation 

Azg geodetic azimuth 

Azg spherical azimuth 

a semimajor axis of earth 

b semiminor axis of earth 

c constant defined in appendix D 

D radial distance from V to T 

E column matrix of residuals 

e eccentricity of earth 

F,L,M,S,A,P,^> functions (Main symbols are used in subscript position to 

indicate the independent variable; for example, a is 
the variable in 

a , n ' 

f flattening of earth 
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horizon uncertainty 

height of surface point above ellipsoid 

arbitrary point in a plane 

normal to ellipsoid 

North Pole 

arbitrary observable 

upper limit of summation 

rectangular matrix of coefficients of observation equations 
row matrix of coefficients of arbitrary observation equation 
radial distance from earth's center (to a point) 
mean radius of earth 

perpendicular distance from earth's polar axis (to a point) 

South Pole 
radar station 

arbitrary point on earth's visible horizon 
column matrix of independent variable 
reduced latitude 

arbitrary location of space vehicle 
e^sin^0 

rectangular cartesian coordinate system (to be defined) 
coordinates of point in X^Y^Z coordinate system 
earth subtense angle 

azimuth of point on earth's visible horizon as seen from vehicle 
elevation of vehicle with respect to spheroid 
elevation of vehicle with respect to sphere 



A small finite change in variable 

8N height of spheroid above geoid 

difference between geocentric and geodetic latitudes 
|e[,|v} column matrices of residuals and corrections, respectively 
e element of E 

components of deflection of the vertical 
0 angle between geocentric radii of vehicle and radar station 

X longitude 

H azimuth of plane of the vertical 

V deflection of the vertical 

p,T unit vectors 

a standard deviation 

0 geodetic latitude 

0' geocentric latitude 

i|f polar angle 

longitude east of origin of datum, Xs - Aq 
S ubscripts: 

H horizon 

1 instrument 

i particular azimuth or point on earth’s visible horizon 

k number of observation equations 

m type of observation 

n particular position of vehicle 

o variable referenced to origin of spheroid 

s radar-station-referenced variables (except where defined differently) 

T total 


k 



sv measTirement from radar station to vehicle 

true true location or value 

V vehicle -referenced variables 

Notations: 

column matrix 
square matrix 
row matrix 


I I 
[] 
I J 


rectangular matrix 
* least -squares estimate 

Matrix exponents: 

T indicates the transpose 

-1 indicates the inverse 

A "bar over a symbol (or group of symbols) indicates a vector. 

An arc over a symbol (or group of symbols) indicates an arc measure. 
A primed function indicates a partial derivative. 


BRIEF REVIEW OF GEODETIC THEORY 


As is well known, the earth's figure is an irregular one. Furthermore, 
the density varies from region to region. Therefore, in order to locate the 
various features on the earth’s solid surface, geodesy relies on two reference 
STirfaces, the geoid and the oblate ellipsoid. The geoid is defined most pre- 
cisely as the equipotential surface of the earth’s gravitation and rotation 
which in the open ocean will coincide on the average with the mean sea level 
and to which the waters on the ocean would tend to conform if allowed to flow 
into very narrow and shallow canals cut throu^ the land (refs. 1 and 2). Thus, 
the geoid coincides, on the average, with the mean sea level in the open seas 
but will rise under the land masses in proportion to the attraction of the land 
mass above mean sea level. 

If the earth’s surface were a perfect ellipsoid, the geoid would almost 
be an exact ellipsoid. However, because of the irregularities in the earth’s 
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shape and the variations in its internal density, the geoid, although smooth, 
is an irregular surface and follows the earth's mean contour. According to 
reference 1 the geoid departs from the spheriodal shape by a few hundred feet 
and in inclination by as much as a minute of arc. (Although there is a dis- 
tinction between the terms spheroid and ellipsoid, the two will be considered 
synonymous in this paper.) 

The vertical axis of a level theodolite at sea level is coincident with 
the direction of gravity at the point of observation. Further, the horizontal 
axis of a level spirit level at sea level is normal to the direction of gravity 
at the point of observation. Therefore, since the direction of gravity at sea 
level is everywhere normal to the geoid, this is the surface to which the geode- 
sist references all of his measurements. As a result, observables such as 
astronomic latitude, longitude, and elevation above sea level are defined in 
terms of the geoid. 

Astronomic latitude and longitude in the conventional sense are defined 
as follows: the astronomical latitude is the angle between the tangent to the 

observer’s plumb line at its intersection with the geoid and the earth’s mean 
equatorial plane and is positive north. The astronomical longitude may be 
defined as the angle between the mean meridian plane of Greenwich (this plane 
contains the normal to the geoid at Greenwich and is parallel to the earth’s 
mean rotational axis) and the plane normal to the earth’s mean equatorial plane 
and containing the observer’s geoidal normal mentioned and is measured positive 
east of Greenwich. To be precise in these definitions, it would be necessary 
to substitute the word "instantaneous" for the word "mean." However, since the 
orientation of the reference ellipsoid relative to the geoid must fulfill the 
condition that its minor axis be parallel to the earth’s mean axis of rotation, 
the earth’s spin axis will be considered to be its mean axis of rotation and 
only the conventional definitions of astronomic latitude and longitude will be 
considered. (See ref. J>.) 

The elevation of a point is defined as its height above the geoid. For 
the purposes of this paper, this definition will suffice. The act of measuring 
an elevation by means of a spirit level assumes that an equipotential surface 
adjacent to the geoid is everywhere equidistant from the geoid. However, this 
is not true and, hence, heights obtained by spirit leveling over great dis- 
tances must be corrected. (See refs. 1, 2, and 4 .) 

An oblate spheroid or ellipsoid of revolution approximates the geoid very 
closely and is the second surface with which the geodesist is concerned. The 
spheroid is a mathematical surface and is by definition smooth. It is, there- 
fore, the ideal surface for geodetic computations and for locating geodetic 
points . 

The geodetic latitude is the angle between the earth’s equatorial plane 
and the normal to the ellipsoid containing the observer’s projected position 
on the geoid. The geodetic longitude is the angle between the meridian plane 
containing the same ellipsoidal normal and the meridian plane of Greenwich 
where both planes are parallel to the earth’s mean axis of rotation. (See 
ref. 3.) 
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If the astronomical latitude and longitude of sin arbitrary point and the 
azimuth of an additional point are desired, it is required only to observe 
these quantities. However, if it is desired to locate a second point relative 
to the first, it is, for the most accurate determination, necessary to connect 
the two points by first-order tri angulation; select an ellipsoid which will 
best fit the region surveyed and define its orientation with respect to the 
first point; then, calculate the geodetic coordinates of the second point on 
the assumed ellipsoid. The resulting geodetic coordinates of the second point 
would, most probably, not agree with those determined astronomically. The 
major reason for this discrepancy is the varying density of the earth from 
point to point on its surface, which would cause the vertical (or direction of 
gravity) at the point to have a direction different from that of the normal to 
the ellipsoid at that point. 

The angle between the vertical and the normal to the ellipsoid is defined 
as the deflection of the vertical. The deflection of the vertical can be 
resolved into two components | and t], which are conventionally positive if 
the inward vertical is directed to the south or west of the normal. (See 
ref. 1.) The component of deviation in the meridian plane q is given by the 
relation 


T] = Astrolatitude - Geodetic latitude (l) 

The component in the prime vertical plane | is given by the relation 

I = (Astrolongitude - Geodetic longitude)cos 0 (2) 

In some texts the definitions of | and q are reversed. 

A second reason for the discrepancies between astronomic and geodetic 
coordinates may be a poorly defined datum or spheroid. The definition of the 
datum of a triangulation network requires eight constants. First, the minor 
axis of the ellipsoid is always defined to be parallel to the earth's spin axis, 
and thereby, the orientation of the triangulation on the ellipsoid can be deter- 
mined through the Laplace equation. This involves first the observation of the 
three constants, the astronomic latitude, longitude, and azimuth and then the 
assigning of values to the two components of the deflection of the vertical 
Iq and tiq. Secondly, values must be assigned to the two ellipsoidal param- 
eters a and f . Finally, the hei^t of the ellipsoid above the geoid SNq 
is defined to complete the definition of a datum. 

The definition of a datum, except for the astronomic latitude, longitude, 
and azimuth, is purely arbitrary. Generally, the geodetic and astronomic coor- 
dinates may be defined as being equal at the origin. However, such an assump- 
tion may not yield the best fit to the geoid. Consequently, it would be 
desirable, in such cases, to define | and q to be different from zero. 

The last three constants |q, Hq, and 6 Nq represent an attempt to 
define the center of the ellipsoid which, ideally, would be at the earth's cen- 
ter of gravity. Because of the irregular shape of the geoid and the uncer- 
tainty in the magnitudes of a and f which best fit the geoid as a whole. 
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this does not hold. Consequently, many independent surveys have been carried 
out which have resulted in several spheroids (European, American, etc.), each 
defined so as to fit best the particular area in question. Attempts to connect 
the various spheroids have shown serious discrepancies. The failure to connect 
the various spheroids is the result of two factors: the true deflection of the 

vertical at the origin for each datum is unknown, and the constants of the best 
fitting world ellipsoid are unknown. 

Stokes has described and developed the theory of a method which, if grav- 
ity observations over all the earth were available , would permit the calcula- 
tion of geold \xndulations at any desired point on the earth's surface. Vening 
Meinesz later extended the theory to show that the components of the deflection 
of the vertical could also be determined from gravity observations. (See 
ref. 5«) When the earth as a whole is considered, the holdings of gravimetric 
data are not by any means complete although enough data exist to make a few 
very good approximations. In the southern hemisphere, observations are sparse 
indeed. In the northern hemisphere, the picture is much brighter but still 
sadly lacking. 

The job of collecting gravimetric data is a gigantic and expensive under- 
taking. Therefore, a more efficient and direct method of determining the best 
fitting world ellipsoid is desirable. 

There are several ways to determine the earth's figure. The most fre- 
quently used approach has been to observe the differences in astronomic lati- 
tude and longitude from place to place and to connect the points of observation 
by triangulation so as to determine the lengths of the degree of latitude 
and/or longitude in different parts of the earth and to deduce the earth's fig- 
ure from the results. Another approach has been through measurement of the 
variation of gravity between the equator and poles of the earth. It was the 
observation of the direction and intensity of the earth's gravitation which led 
to the discovery of irregularities in the earth's figure and internal structure. 

There have been a number of attempts to derive an ellipsoid which will 
best fit (as a whole) the earth's figure. Until 1959 j the international ellip- 
soid (a = 6 378 388 meters; f = l/ 297 ) was generally accepted as the truest 
representation of the earth's figure. However, in reference 6, observations 
of the Vanguard satellite were shown to indicate a flattening of 1/298.3- 
These findings resulted in a rash of additional derivations of the earth's 
ellipsoidal parameters. However, there is still disagreement. 

Once the parameters of the best fitting ellipsoid have been established 
it would be necessary either to change the spheroids of the various geodetic 
systems and recompute the station coordinates or to derive formulas for con- 
verting from the old datum to the new. The two alternatives differ only in 
the extent to which the computations are carried. For the first alternative 
the coordinates of all geodetically established points would be changed. For 
the second, a convenient formula would be provided for converting the coor- 
dinates of any desired point on the old spheroid to those cf the new. 
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If it is desired to change a spheroid, one approach would he first to 
redefine a, f, |o^ tJq, and SNq and then compute the coordinates of the 
observed station of the triangulation network again. There are available sev- 
eral methods, of which the methods of Vening Heine sz and de Graaff-Hxmter are 
two that allow for a change in the datum without going throu^ the rigoroxis 
confutation of the triangulation just outlined (refs. 1 and 7)* Essentially, 
if At)q, A|q, A6Nq> Aa, and Af are the changes to be added to tJq, |q, 

&Nq, a, and f, respectively, the corresponding change in the geodetic lati- 
tude A0s^ longitude AAs, and station height Z!hs may be obtained directly. 
The de Graaff-Hunter method was chosen for this study. In functional notation 

A0S = L(Aa,Af,A|o,ATiQ,ASNor 

/SKs = M(Aa,Af,A|o,ATio,A8No) > (3) 

= S(Aa,Af,A|o,ATio,A6No)_^ 

Assming the availability of all desirable data, to establish a truly 
world geodetic system and relate all existing datums to it by the previously 
discussed procedures is an overwhelming task. Therefore, a method which would 
allow the dat\am to be adjusted in addition to deriving a truly earth-centered, 
best-fitting world ellipsoid would be of great importance. The development of 
such a method is the purpose of this paper. Because of the magnitude of the 
problem, however, the present report proposes to develop the concept according 
to certain assumptions and apply it to a particular case. Some of the asso- 
ciated problems and areas for addltonal work are pointed out as well as indi- 
cations of the merit of the method. 


ASSUMPTIONS AND GENERAL THEORY 


If an artificial earth satellite is assumed to be at a point Vn in space 
and Rv,n units from the earth’s center at time to, as is shown in figure 1, 
and, by some suitable optical means , a point T^, on the earth's visible hori- 
zon, is observed, the vector V^^Tj^ (assuming no atmosphere) is the line of 
sight and is tangent to the earth at Ti. If the earth's figure is assumed to 
be a smooth oblate ellipsoid, the vector is perpendicular to the n ormal 

Ni to the ellipsoid at Tj^. The angle an,i included by vectors VnTi and 
in the plane of observation is defined as the Ttli subtense angle of the 

earth. 

It is obvious that, in addition to being dependent on the size a nd sh ape 
of the earth, otn,i will vary in magnitude as the plane containing VnTi and 

Rv,n rotates about Ry,n ^<i S'S the location of Vn changes with respect to 

the earth's center. Further, if the earth is assumed to be a smooth oblate 
ellipsoid of revolution, due to the symmetry of the ellipsoid, On,i is 
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Figure 1.- Optical observation of a point on earth* s visible horizon 

from a satellite. 

independent of the vehicle* s longitude if 3^ is the azimuth of 

the plane formed hy V^Ti and Ry^n measured in a plane normal to Ry^n 
clockwise from the meridian plane of Vn and if the spherical coordinates of 
Vn are Ry^n^ ^v,n-» ^,n^ where 0y^n is the geocentric latitude of 

the earth subtense angle for the i^^ point of tangency and the n^^ position in 
space is written as 

^^i “ ^(^^^>^v^n^0v^n^Pi) 

In this analysis^ the assumption is made that the lines of sight are truly 
tangent to the solid surface of the earth and are not intercepted by clouds^ 
haze^ or atmospheric scattering. Effects of atmospheric refraction are 
neglected. 

The earth's solid surface is an irregular figure. However^ if it is 
assumed that the elevations of various features on the earth's surface vary in 
a random manner about some mean elevation, the summation of all the elevations 
about that mean over any path would be zero. That is, if the earth's mean 
shape is assumed to be that of an oblate ellipsoid of revolution, there is an 
ellipsoid with semimajor axis a and flattening f. If this ellipsoid is 
placed so that its geometric center is coincident with that of the earth (the 
earth's geometric center and center of gravity will be taken to be synonymous), 
the s\m of all the elevations about it will be zero. Thus, if a large number 
of subtense angles are observed from an artificial satellite at different 
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positions in space and the geocentric radii and latitudes are observed from 
some earth-based observation station, the datum to which the observation sta- 
tion is referenced could be Improved by the method of least squares. 

An ellipsoid determined in the manner described here will represent the 
earth’s mean geometrical shape. Further, the resulting fit to the mean geo- 
metrical shape will be a function of both the number of observations and the 
particular surface features observed. For example, if observations are made 
predominately in a mountainous region, the resulting ellipsoid would be larger 
than desired. On the other hand, however, if observations are made only to the 
surface of the oceans, the resulting ellipsoid would be smaller than desired. 
(See previous section.) Thus, if observations should be made to all regions 
of the earth’s surface, the ratio of the number of observations on a particular 
feature, such as mountain ranges and oceans, to the total can be assumed to be 
proportional to the ratio of the area of the feature to the total surface area 
of the earth, amd the resulting ellipsoid would represent a best fit to the 
earth’s mean geometric shape. Axi error exists in using horizon measurements 
in mountainous regions, because of the fact that the lines of sight from a 
satellite to the horizon may be \xnable to reach the ground in valleys. Because 
the oceans and seas cover approximately 2/3 of the earth’s surface and the 
mountainous regions comprise a very small portion of the total surface area, 
this error would probably influence the resulting ellipsoid only slightly. 

The two spherical coordinates of the vehicle Ry,n 0v,n equa- 

tion ( 4 ) are dependent on the geodetic coordinates of the observation point and 
the observations necessary to fix the location of the vehicle in space. In 
functional notation 

0v,n ~ ^ . . 

^v,n “ ^(^j^^0s#^s^^s^Om,n^*^m-l,n^ • * 

where the values of m (m=l,2,5^***) represent the particular observa- 
tions necessary for the purpose of fixing the position of the vehicle in space. 

As explained in the section on geodetic theory, the geodetic coordinates 
of the observation station in equation (5) are not independent of the origin 
of the datum. Thus, since the datum is arbitrarily defined so as to obtain a 
best fit of the ellipsoid to the surveyed area, the spherical coordinates of 
the vehicle must be considered to be in error. The error in each is given by 
its total differential. Consequently, the error equations are 


^°l,n 
^Ol,n)) 
■^Ol,n) 


( 5 ) 
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d0'=Mda + |idf + |l 


d0s + dAg + ^ dhs + y — dOm n 
SXs ^ 5hs ® L SOjn n ^ 

m=l ’ 




> ( 6 ) 


m=l 

P 


■ S ^ " i “ * I; * I; ^ * Z ^ 

m=l 

Those variables Ojjj jj which are to be observed by some suitable method 
may change with the ellipsoidal parameters a and f and the geodetic coordi- 
nates of the observation station. Thus, if it is assumed that the error in each 
dOm n given by the relation 


dO, 


So, 


m 


So, 


So, 


So„ 


So„ 


Sa 


^ da + df + d0g + dAs + dhs (7) 


Sf 


S0 


SAf 


She 


upon substituting equation ( 7 ) into equation (6) and collecting coefficients of 
like terms, error equations for the spherical coordinates of the vehicle as 
functions of the errors in the ellipsoidal parameters and geodetic coordinates 
and height of the observation station are 


d0* = O* da 

"v,n a,n 

<i\f,n ~ ■^,n da 
d^,n “ ^a,n 


+ ^>f^n df + d0s + 

+ Af^n df + A^g^n d0s + A^^^^ 

+ df + d0s + PXg,n 




+ %,n 


dAg + ^hg^n ^s\ 


P^3,n '“s 


(8) 


where the A*, and P* coefficients represent the sums of the partial 

derivatives of the particular function with respect to the variable indexed. 


The differential operator d of equation (8) is equivalent to the A in 
the relations of equation ( 3 ). The mathematical difference in the two, for 
example, is that da represents an infinitesimal error in the variable whereas 
Aa represents a small but finite error in the variable. However, it is assiimed 
that the coefficients of the unknowns are linear over a finite range of the 
variables. Thus, the two mathematical representations are taken to be equal, 
and, therefore, the error equations for the spherical coordinates of the vehicle 
can be determined as functions of a change in the datura. 
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The relations of equation (5) are linear functions of their variables. 
Therefore , 

A0S = ^a^ ^f^ ^^'Ho ^^^o ^ 5 N^^o"^ 

Ms = <Aa + + M^AtIo + M|A|o + M^igASNo > (9) 

Mg = S^Aa + + S;|At1o + sJa^^ + S^6No_^ 

where primed-fanction names represent the coefficient of the corresponding 
unknown^ which is given as a subscript. Then^ after substituting equation (9) 
into equation (8)^ 


^;,n = ('^;,n + + (^f,n + 4*03, n 44s,n 

+ 4*hg,n)^ ^ (4*0s,n + ^ 44s,n)^^o + (4*iZSs.n 

^i*hg,n)^o + (4N*0s,n 4N*As,n 4N*hg.n)^^^o 

= (4,n ^a.^s,n ■•■ ®aAhs,ny^ (Af,n ■*■ ^f^g,n ■*■ ^%^Ag,n 

■*■ ®f^s,n)^ + 4 '^hsjn)^o + ^|^As,n 

+ (^6N4s,n ^ 5 N^As,n ^6lAas,n]^^o 

AEV,n = (P;,n + H- 

sXg,n)^ (44s, n + 44s,n + 44g,n)Ario + (44s, n ^ 44s,n 

+ 44 s,n)A^o + (^^4g,n + 4 N 4 s,n + A 5 No 


> ( 10 ) 


If the coefficients of equation (8) can be determined and are linear in 
the corresponding variables^ the change in the spherical coordinates of the 
vehicle due to a change in the datiim can be computed by equation (lO). Thus^ 
only the error equation for remains to be found. 
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Since the fiinction F of equation {k) may not be linear in its variables, 
it is linearized by expanding it into a Taylor series and ignoring all terms 
containing powers of the variables higher than one. Thus, if 


dpi = 0 

for all the remaining variables of F 


°ha,i - (^)n,i 



da + 



+ 


^F \ 

,^^v,n/n,i 


dRy,n 


Sf 




d0. 


v,n 


v,n/n,i 


(11) 


where a partial derivative of a f\mction enclosed by parentheses and subscripted 
indicates that it is to be evaluated for a particular azimuth at a point . 

Equation (u) can be interpreted as stating that 

Aan,i = Error in - (F)n^l = 0 ^^! (true) - o^^i (approximate) 



Substituting the necessary expressions from equation (lO) into equation (12) 
and collecting the coefficients of like terms gives the following error equa- 
tion for 


^ + W)n,l " * (''^^,1 {4o)n,l 


d5W. 


(15) 


where F* is the coefficient for the variable subscripted. 
Generali zing further , let 

= ^,i 

N “ «)n,i' 

|u| = da^ df^ ^To^ ^^o-^ d6NQ 






where 


( 14 ) 


k = 1, 2, . . . 
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is a particular observation. Then, in matrix notation 


€k - Lqkjlul (15) 

Equation (15) represents only one ohse 2 rvation. For a least-squares deter- 
mination of the unknowns to he valid, it is necessary to have a large number of 
observations. In the theory presented so far, two varying indices i and n 
have been used. Thus, if |e| is chosen to represent a column matrix of k 
elements and ||q| 1 is chosen to represent a k by j rectangular matrix of the 
coefficients, where j is the number of unknowns and v is a colvimn matrix of 
corrections, equation (15) can be rewritten to read 

INIM = 1®! + M (16) 

which after requiring the sum of the squares of v to be a minimum becomes 

= iioiriE} (17) 

The product ||Q|P[|Q|| is a square matrix and can be inverted. Therefore, 

iu*i = [[ii«ii''ii«ii] (18) 

is the matrix representation of the least-squares solution for the unknowns 
(ref. 8). A similar solution would be obtained for weighted observations. 

As noted previously, atmosphere refraction has been neglected in develop- 
ment of the theory. If equation (4) is referred to, it will be seen that the 
observed subtense angle can easily be corrected for atmospheric refraction pro- 
vided a suitable correction is available. 


THE EXPRESSIONS FOR a AND ITS ERROR EQUATION 


The theory given in the previous section was very general in that func- 
tional notation alone was used. In this section the expressions for On^i 
and its error equation are given. 

The Expression for a 

If a satellite is located at a point Vn in space and if by some s\ii table 
optical means (assuming no atmospheric refraction) a point T^ on the visible 
horizon is observed, the subtense angle On,! is given by the functional rela- 
tion of equation (4). Fiirthermore , if the earth is assumed to be a smooth 
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oblate ellipsoid, the line of si^t V^Ti is tangent to the earth at Ti and 
perpendicular to the normal Nj^. Then, it is shown in appendix A that the 
desired expression for is 



The positive sign is chosen so that is a minimum when Pn,i zero. 


The Error Equation for a 

Before il>. A, and P can be derived, it is necessary to decide what 
method will be used to fix the position of the satellite. Because of the high 
accuracy obtainable, radar was chosen for this paper. Further, it is assumed 
that perfect measurements are made. 
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There are several techniques for fixing the instantaneous position of a 
space vehicle by radar observations in use today. Since it is impractical to 
discuss each technique and its individual application to the theory of this 
work, the conventional measurements of range azimuth AZg^j^, and the 

elevation above the horizon are taken to be observables. 

The azimuth and elevation are geodetic azimuth and geodetic elevation. 

Thus, if the datum should be changed, Azg^j^ and 7 ^ must necessarily change. 

On the other hand, however, the range Rgy^n ^ pure measurement which 

is solely a function of the radar accuracy. As a result, derivatives of Rsv,n 
are taken to be zero. 

With the assumptions stated herein in mind, it is shown in appendix B that 
the coefficients of the expressions in equation (8) can be expressed as follows: 
For d0;^^ 


<D 


a,n 


f,n 


^^v,n ^s ^^s ^ ^v,n ^®n ^v,n ,n 




8Ng 8a 


80 n 8a 


8Az 


s ,n 


8a 


^v,n/^s ^s ^s\ ^^v,n ^n ^^v,n ^^s,n 


80’ \8Ns 8f 8f 


80 n 8f 


8az, 


8f 




_ ^v,n ^9n . ^v,n ^^s,n 


80n 8As 8Az 


s ,n 


8Ac 


hg,n 


.. ^v.n ^s , K,n ^Qn ^ ^v.n ^^s.n 


80' 8hg 


80 n 8h^ 8 Azc: 


8hc 


E>* = ^^s ^ ^s\ ^v,n ^^n ,n ^^s,n i 

0 s. n - ^ 0 . 80s 80s/ 80 n 80n 8 Azs^n ^ 0 s [ 




( 21 ) 
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For dA, 


For 


v,n 




^SY,n ^ ^sv,n ^sv,n 

hQj^ be. 




bd' 

*^v,n 


S jU 


^,n--»An a0;^; 


^^sv,n ^^sv,n ^®n ^^sv,n ^^s,n 

4 * - ' " - 4 “ — 


SlAZg^jj Sf 


j^y - $* ,n ^‘^sv,n ^0n ^ ^^sv,n ,n 

%,n = 0s-”l0^ a«n ^ aAzg^n 30, 


Ai„n = 1 + < 


(VvA ^9 SiA.z„ _ 

sv,n ^ sv,n n sv,n s,n 


^s>n SAzg^n 


ckAAg 

is.n 


a’ = o' ^^sv ,n . ^^sv,n ^^sv,n ,n 

hs,n h.,n ^^s,n ^ 




P' = 0 * ^liH + A* ^^r,n ^^v,n ^^s,n 

^7k ^ ^-s,n ^ 


(22) 




P' =$' 5 i£I + a' K , ^^s,n 

f,n f,n f,n ^Zg^n 


= o’ ^ ^ ^Ry>n ^^s,n 

0s.n ^s,n a0g ^Zs,n ^s 


> ( 23 ) 


p ' = oJ + (aJ _ i\ JVjiL. + M 

As.n As.n ^ As,n BAs ^Zg^„ SAg 


= o' t ^^,n ^7n ^^s,n 

^;,n ^sv,n ^^s,n 


tig;ii h„ ,n 
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Firrther, the coefficients of equation (l2) are shown to he given by 


. ^,i ^^i,i ^ ^,i ^B„,i ^ ^,i ^n,i 
^'n^i ^n,i ^ ^®n,i ^n,i 


(M.\ = ^>1 ^,i ^,1 ^^,i ^ ^,i ^n,i 

wf/n,i ^n,i ^®n,i ^n^i 





^,i ^,i ^,i ^\,i ^ ^,i ^n,i 

^n,i ^v,n ^®n,i ^v,n ^n,i ^v,n 



n 


A 


^n.i , ^,1 ^^n,i , ^,i ^n,l 

^^,i ^®n,i ^^,n ^n,i ^^,n^ 


(24) 


The error equation for Oji i given in a general form by equation (ij). 

As stated previously^ the coefficients of equation (15) are functions of those 
of equations (8)^ (9), and (l2). 

For brevity^ in this section the coefficients of equations (8) and (l2) 
were given in equations (2l)^ (22)^ (23)^ and (24) as functions of indicated 
partial derivatives. The indicated operations have been executed and the 
results are given in appendix C. To obtain the algebraic expressions for equa- 
tions (21)^ (22)^ (23)^ and (24)^ it is required only to make necessary 
substitutions . 


To complete the list of expressions for the solution of equation (l3)^ it 
is only necessary to rearrange the equations due to de Graaff -Hunter to the 
form given by equation (9)» This has been done and the results are given in 
appendix D. 

The foregoing has shown that the error equation for ^ is determined 

in terms of quantities which are known^ are assumed to be observed^ or can be 
approximated. Thus^ it is mathematically possible to derive by the method of 
least squares an earth-centered geodetic datum from earth subtense angles 
observed from an orbiting artificial earth satellite. It is now desirable to 
mention ways in which the theory can be applied and some of the conditions 
necessary for obtaining valid results. The former is the subject of the next 
section and the latter will be discussed in a subsequent section. 
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APPLICATION 


There are several ways in which the theory of this paper can he applied, 
hut each will have one point in common with the other. All methods will involve 
the direct observation of some angle either between two points on the earth's 
horizon, between a reference line and points on the horizon, between a refer- 
ence star and points on the horizon, or between some other reference and points 
on the horizon. The approach chosen for a particulsir application will be 
largely dependent on the type of instruments available and their acc\iracies. 

For comparison and discussion this section will consider two applications. 
One is simply that of observing the subtense angle directly and another is that 
of observing angles from a polar star to points on the earth's visible horizon. 

There is no need to discuss the first approach in any greater detail. The 
assvuBption of being able to measure the earth's subtense angle by some means 
either directly or indirectly was the basis for the derivations which precede 
this section. The only drawback to such an approach is that the geocentric 
radius is not a physical reality, but mathematically defined. Therefore, 
Instnamentation would be required to operate in such a manner as to allow the 
precise determination of the direction of the geocentric radius. 

The approach of observing directly the angle between the direction of a 
polar star and the line of sl^t to a point on the earth's horizon is less com- 
plicated from the standpoint of instrumentation than that of observing subtense 
angles directly. Since the direction of a polar star would be known, it would 
only be reqiiired to observe the direction of a point on the horizon relative to 
the star. To utilize the "polar angle" determined in this manner, only slight 
modification of the procedure developed in this paper is necessary. 


Consider figure 2 in which the Z'-ajcls is considered to pass throu^ the 
vehicle's center and is positive in the direction of a polar star. If the star 
is considered to lie directly above the earth's north pole, the Z'-axis is 
always parallel to the earth's spin axis regardless of the location of the 
vehicle, since the star is essentially infinitely distant. As a result, Z' 
would translate about the earth's spin axis and, if +n,i polar angle. 


0v,n °^,i " ^v,n “n,l cos ( 25 ) 


Then, 


^+n,i = ^ <,n + ^ ao^,i 
^v,n ^,i 


(26) 
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Figure 2.- Geometric definition of polar angle and its relation to suttense 
angle, azimuth of observation, and spherical latitude of a satellite. 

where 

_ cos cos + sin sin cos 

^v,n 'fn^i 

> ( 2 ?) 

^n,i ^ _ sin sin + cos cos cos p^ 

aon^i “ sin 

Therefore , since the expressions d0y and <ian,i known and the partials 

of equation (25) are given hy equation {2'j) , after substitution the procedure 
becomes the same as that for the subtense angle. 

In the absence of real data a computer program was written for the purpose 
of simulating data. The results were then used to determine the resulting 
earth-centered datum. 

The method of simulating measurements was briefly this: A datum was 

defined at a selected latitude and longitude along with the coordinates of a 
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radar station; various positions of the vehicle were assumed in terms of radar 
measurements, and the approximate subtense angles were computed for certain 
values of 3j_; errors were then assumed to exist in the datum, and the true 
coordinates of the radar station were determined by the method of 
de Graaff -Hunter (appendix D); finally, the assumed azimuth and elevation of 
the vehicle were corrected by equations (Bl4) and (BI 5 ) in appendix B, and the 
true subtense angles were found for the assumed values of 3j: . A similar pro- 
cedure was used for the polar method. 

The coordinates of the origin of the datum and the radar station and the 
errors in the datum are given in table I. For simplicity, the vehicle was 
assumed to he moving away from the earth at a constant elevation of approxi- 
mately azimuth of 225°. Observations were assumed to be ini- 

tiated when 


Pi = 0.0° 


and 


Rsv = 200 000 meters 

where Pi was increased in steps according to the relation 


Pi+1 - 100 


whereas Rgy increased in increments of 200 000 meters. One observation 

was assumed to be taken for every new value of 3j_ and Rsy* table II the 

resulting angles and residuals are given for the assumed true and approximate 
subtense and polar methods for each geocentric latitude and longitude. Fifty 
observations were taken and a solution was found after each, beginning with the 
fifth observation. 


TABLE I.- DEFINED CONSTANTS FOR TEST CASE 


Variable 

0, radians 

A; radians 

h, meters 

a, meters 

f 

^ ^ radians 

T), radians 

Origin of datum 

0.52359870 

-1.3089967 

10.0 

6378388.0 

336.70834x10-5 

-9.692736x10-6 

1.4515551x10-5 

Radar station 

.58067851 

-1.3835296 

.0 



2.4242184 

2.0000000 

Change in datum 



- 8.0 

-223.0 

-1.4673498 

4.8484578 

-1.2120342 
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TABLE II.- TRUE AND APPROXIMATE SUBTENSE AND POLAR ANGLES AND THE CORRESPONDING RESIDUAI^ FOR THE TEST CASE 


0 v,n^ radian 


radians 

radians 

Approximate 

True 

^n,i 

Approximate 

True 

^n,i 

0.5659^264 

6 541.37880 

1.33965063 

1.33965151 

0.08770755 X 10-5 

0.79708833 

0.79711264 

0.24315266 X 10 -^ 

.55480848 

6 712.47816 

1.24795322 

1.24795581 

.25974485 

.87816855 

.87819012 

.21578650 

.54416163 

6 885.15547 

1.17941563 

1.17941847 

.28401545 

.93747975 

.95750010 

.20355788 

.55397802 

7 059.25635 

1.12312724 

1.12512990 

.26640128 

.98578545 

.98580305 

.19593455 

.52423424 

7 234.67654 

1.07481548 

1.07481781 

.23209582 

1.02724484 

1.02726388 

. 19040745 

.51490762 

7 411.36098 

1.03226408 

1.03226598 

.19020476 

1.06399216 

1.06401077 

.18605854 

.50597644 

7 589.20282 

.99414064 

.99414209 

.14475538 

1.09728131 

1.09729955 

. 18246851 

. 49741994 

7 768.12245 

.95956622 

.95956720 

.09775206 

1.12792277 

1.12794071 

.17941519 

.43921835 

7 948.04716 

.92792249 

.92792299 

.05027029 

1.15647264 

1.15649051 

. 17676648 

. 48135288 

8 128.91018 

.89875341 

.89875344 

.00291677 

1.18532969 

1.18334714 

.17443574 

. 47380581 

8 310.65023 

.87171028 

.87170984 

-.04595416 

1.20878971 

1.20880695 

.17236107 

. 46656035 

8 493.21103 

.84651871 

.84651781 

-.09013019 

1.25307787 

1.25509492 

.17049487 

. 45960066 

8 676 . 5'^85 

. 8 ?? 9578 i 

.82295645 

-.13548122 

1.25636925 

1.25658612 

.16879836 

. 45291182 

8 860.59189 

.80084630 

.80084450 

-.17992547 

1.27880223 

1.27881896 

. 16723843 

. 44647975 

9 045.52015 

.78005308 

.78003085 

-.22340959 

1.30048792 

1.30050449 

.16578592 

. 44029124 

9 230,68488 

.76039049 

.76038783 

-.26589677 

1.32151635 

1.32155280 

,16441452 

.43433382 

9 416.64851 

.74180942 

.74180634 

-.30735944 

1.54196140 

1.54197772 

. 16310018 

.42859579 

9 605.17636 

.72419570 

.72419222 

-.34777488 

1.36188406 

1.36190024 

. 16182072 

. 42306611 

9 790.23612 

.70746755 

.70746549 

-.58712255 

1.38155509 

1.58135116 

. 16055561 

.41773441 

9 977.79799 

.69155248 

.6915^23 

-.42538246 

1.40035701 

1.40057294 

. 15928581 

.41259096 

10 165.83383 

.67638758 

.67638295 

-.1+6253444 

1.41898555 

1.41900154 

.15799372 

.40762655 

10 354.51826 

.66191624 

.66191126 

-.49855760 

1.43725084 

1.43726650 

. 15666509 

.40283257 

10 543.22684 

.64808811 

.64808278 

-.53543026 

1.45517844 

1.45519597 

.15527905 

.39820085 

10 732.53759 

.65485798 

.65485231 

-.56715004 

1.47279003 

1.47280541 

.15582805 

.39372376 

10 922.22905 

.62218514 

.62217914 

-.59963401 

1.49010599 

1.49011922 

.15229783 

.38939407 

11 112.28216 

.61005275 

.61002643 

-.63091897 

1.50713594 

1.50715101 

.15067755 

.38520498 

11 302.67847 

.59836751 

; .59856070 

-.66096170 

1.523899:4 

1.52591402 

.14895761 

.58115006 

11 495.40105 

.58715840 

.58715151 

-.68973931 

:. 54040478 

1.54041949 

.14712977 

.37722525 

11 684.43394 

.57637817 

i .57657100 

-.71722952 

1.55666232 

1.55667685 

. 14518707 

.57541885 

11 875.76199 

.56600112 

1 .56599369 

-.74341100 

1.57267968 

1.57269400 

.14512386 

.36973146 

12 067.37125 

.55600384 

.55599616 

-.76826562 

1.58846346 

1.58847755 

.14095579 

.36615597 

12 259.24861 

.54656475 

.54635683 

-.79176875 

1.60401905 

1.60405290 

.13861977 

.36268757 

12 451.38156 i 

.55706596 

.53705585 

-.81390952 

1.61935084 

1,61956446 

,13617394 

.35932167 

12 643.75854 1 

. 52808310 

.52807475 

-.83467105 

1.63446228 

1.63447365 

.13359768 

.35605398 

12 836.36856 i 

.519^512 

.51959657 

-.85404059 

1,64935607 

1.64936915 

.13089154 

. 55288040 

15 029.20127 ' 

.51101422 

.51100550 i 

-.87200791 

1.66403410 

1.66404691 

.12805722 

.34979706 

13 222.24689 

. 50289571 

.50288682 

-.88856524 

1.67849769 

1.67851020 

.12509756 

.34680028 

13 415.49623 

.49505592 

. 49502688 

-.90570756 

1.69274753 

1.69275974 

. 12201642 

.34388659 

13 608.94072 

.48742210 

.48741292 

-.91745265 

1.70678383 

1.70679572 

.11881872 

.34105268 

15 802.57201 

.48004235 

.48003306 

-.92974:26 

1.72060630 

1.72061785 

. 11551031 

.33829541 

13 996.38247 

.47288556 

. 47267616 

-.94063716 

1,73421425 

1.73422544 

.11209795 

.53561179 

14 190.36460 

.46594132 

.46593182 

-.95012715 

1.74760655 

1.74761741 

. 10858925 

.33299901 

14 384.51171 

. 45919986 

.45919028 

-.95822120 

1.76078182 

1.76079232 

.10499261 

.33045435 

14 578.81689 

.45265204 : 

.45264240 

-.96493241 

1.77373831 

1.77574844 

.10131710 

.52797527 

14 773.27418 

.44628929 

.44627958 

-.97027703 

1.78647399 

1.78648375 

.09757244 

.32555931 

14 967.87751 

. 44010350 

. 44009376 

- . 97427444 

1.79898660 

1.79899597 

.09376891 

.32320416 

15 162.62126 

.434087U 

.43407734 

-.97694713 

1.81127362 

1.81128262 

.08991726 

.32090760 

15 357.50020 

.42823295 

. 42822516 

-.97832064 

1.82333237 

1.8235^97 

.08602865 

. 31866752 

15 552.50907 

.42253428 

.42252450 

-.97842351 

1.83515994 

1.83516817 

.08211453 

.31648193 

15 747.64323 

. 41698475 

. 41697498 

-.97728717 

1.84675333 

1.84676115 

.07818659 


RESULTS AND DISCUSSION 


In figure 5 the five arbitrarily defined constants of the datum are plotted 
as a function of the number of observations. The results for both the subtense- 
angle case and the polar-angle case are shown. 

For both cases the results converge to the true values quickly^ as is 
expected. Theoretically^ since no random errors have been assumed to exists 
convergence should be instantaneous. However^ because of certain approxima- 
tions and computer noise ^ instantaneous convergence will not result. Of the 
two methods^ the polar-angle method appears to give the better results. 
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One peculiarity is noticed in the case of the subtense angle. A sharp 
spike occurs in these curves at the sixth observation and is most pronounced 
in the A|q curve. This oddity is not the result of computer noise. 

There are two possible reasons for the spike effect. One is a poorly con- 
ditioned matrix, and the other is a possible breakdown in the equations at the 
sixth observation for this example of the subtense case. 

A detailed study to determine which possibility, if either, is reasonable 
is beyond the scope of this paper. However, two cases were run to determine 
the effect on the spike of increasing the number of instantaneous observations 
at a given point. One case assumed 20 instantaneous observations, and the other 
assumed 200 instantaneous observations at each point. The spike disappeared in 
the latter case and was still present in the former. Since no discrepancies in 
the remaining determination are noted, the spike effect will not be discussed 
further . 


It will be recalled that in the derivations which preceded this section, 
60s was considered to be a variable. In reality the maximum value of 60s is 
approximately 12' of arc at 45° latitude and should remain relatively constant 

for a small change in latitude ^that is, ^ ^ when A0g is small j. 

To determine the effect of 60g on a determination, two cases were run in 
which 60g was first assumed to be a constant by eq.uating its partials to zero 
and then a variable. Further, the approximations 


sin 60g « 80£ 


cos 60s 1 

were assumed valid in both cases. The results are shown in table III for the 
first determination of each case only, since it is representative of the 
remaining determinations. 


TABLE III . - CALCULATED CHANGE IN DATUM FOR TWO CONDITIONS OF 


Assumption 

meters 

— 

Af 

radians 

At|q, radians 

A 6 Ns> meters 

60 s = Const Sint 




h 

4.5626973 

Sin 60 s 60 s 

Cos 60 s 1 

60 s = Variable 

-71.126503 

-0.10408976 X 10"^ 

5,4794727 X 10"5 

-0,14433988 X 10"^ 

Sin 60 s « 60 s 
Cos 60 s ~ 1 

-221.551^3 

-.14628869 

. 51204551 

-.12143458 

-7,8855827 
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When the results of table III are compared with the desired results given 
in table it is seen that has a large influence on the results, which 

was not expected, since 50g is so small and increases slowly. When 60s 
considered to be constant the results were unacceptable although systematic. 
When approximations for the sine and cosine of 60s were substituted for their 
true values, the results were considerably better although not as accurate as 
hoped for. As a result, it is concluded that for best results, 60s cannot 
be assumed constant, and it is necessary to use the true sine and cosine of 
60s in the methods of this paper. 

Since it has been intended to show that the method is mathematically fea- 
sible and to show some of the conditions under which best results are obtained, 
results discussed thus far were computed in double precision on an electric 
data processing machine. Although it is not the purpose of this paper to do a 
statistical analysis of a practical application of the method, it was decided 
to generate some random numbers which could be added to the observed angle for 
the polar method. Thus, it could be determined if the method will converge or 
not with a large standard deviation of measurement. 

The variance of the subtense angle is given by the relation 


CTip — 0j + 0JJ 

where is the variance of the instrument and the variance of the horizon 

uncertainty is given hy the relation 






- R, 




(from ref. 9)- For the case of the polar angle^ the relation for Ojj differs 


from that for the subtense angle in that cr^ must he multiplied hy the expres- 
sion for ^n,i/^*^,i (see eq. ( 27 )) to satisfy the condition that he a 

minimum for Pi = § radians and for 


Pi = ^ n 


radians. After a relation for 


the variance of the polar angle has been determined, random numbers can he 
generated for the purpose of testing the convergence of the method. 


The case tested was identical to the one used throughout this paper with 
two exceptions. First, the incremental increase in Rgy was reduced to 

2 kilometers to keep Rgy small. Secondly, the azimuth A-Zg^n eleva- 
tion were increased in steps of -10"^ radians. The incremental changes 

for each variable were applied simultaneously, and one observation was assumed 
for each point in space. 
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The test case was computed for two instrument variances. In hoth cases ^ 
the horizon uncertainty H was assumed to be 1.6 kilometers. For the first 
case, Oj was taken to be 5 X 10"5 radians, and for the second case, crj was 
5 X 10**^ radians, (in ref. 10, a horizon scanner is described with an accuracy 
of 0 . 1*^ of arc . ) 

The results in both cases were encouraging in that the sigmas for each 
unknown do, in fact, converge. However, convergence was slow. The results for 
both cases in which 50 000 observations were taken are shown in table IV. It 
should be noticed that the standard deviations for Af, ^^o^ ASNq 

are still larger than the desired values (table l). 


TABLE IV.- LEAST-SQUARES ESTIMATE AND STANDARD DEVIATION FOR 50 000 OBSERVATIONS 


standard deviation 
of 

instrument error, 
radians 

Result 

Aa, meters 

Af 

A|q, radians 

radians 

A6 Nq, meters 

5 X 10-5 

Estimate 

-0286.2080U 

-0.22545918x10"^ 

0.06674525x10"^ 

-0.05669926x10-^ 

0190.94084 

a 

0069.21024 

0.09471759x10"^ 

o.ii 452682 xio "5 

0.16108228x10-^ 

0124.95359 

5 X 10 '^ 

Estimate 

-0512.92549 

0.05580511x10"^ 

-0.26521512x10'^ 

0.08405173x10"^ 

0069.04514 

0 

0588.95248 

0.47578258x10"^ 1 

0.45795054x10"^ 

0 . 51414847 x 10 "^ 

0548.57976 


The results discussed here do not by any means constitute a complete error 
analysis of the method. No effort has been made to incorporate observations 
from more than one radar station. Further, it is not known how multiple horizon 
observations at each space point will affect the rate of convergence. The 
results do, however, demonstrate that the answers will converge for large var- 
iances in the observations when sufficient data have been collected. 


CONCLUDING REMARKS 


A theoretical method which can be used to derive an earth-centered goedetic 
datum has been derived and tested by mathematical simulation. Although a prac- 
tical application of the method must necessarily include instrument errors and 
corrections for refraction and aberration, in this simulation such errors and 
corrections were omitted. 

The results of the simulation for a test case proved it to be mathemati- 
cally feasible. Further, when large variances in the observations were assumed, 
the method was found to converge, although slowly. Even so, additional statis- 
tical studies of the method are desirable. As for immediate usefulness, some 
of the req^uired relations may have application in other areas of research - 
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such as, accuracy penalties resulting from assuming the earth to he spherical 
rather than oblate in navigation schemes and horizon - uncertainty studies. 

Tests conducted to determine if the difference between the geodetic and 
geocentric latitudes 60g could be considered to be constant for small changes 

in latitude and sufficiently small to allow its sine and cosine to be approxi- 
mately 50g and 1, respectively, gave unacceptable results. VThen the trigo- 
nometric functions were approximated only, the results were acceptable although 
the accuracy was not as good as hoped for. It is therefore concluded that for 
best results, no approximations can be made in 60g, but acceptable results can 
be obtained by approximating the trigonometric functions of 60g only. Finally, 
under no conditions, could 80g be considered constant. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., December 10, I 965 . 
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DERIVATION OF a 


As explained in the text, if a satellite is assumed to be located at a 
point Vn in space and if by some suitable optical means (assuming no atmos- 
pheric refraction) a point on the visible horizon is observed, the subtense 

angle is given by the functional relation of equation (!<-). Furthermore, 

if th e earth is assumed to be a smooth oblate ellipsoid, the line of sight 
VnTi is tangent to the earth at Ti and perpendicular to the normal Nj^. 

These two assiunptions are the basis for the derivations which follow. In addi- 
tion, the vehicle is considered to be frozen in space with respect to the earth, 
since is independent of longitude as will be shown subsequently. 

Consider figure A-1 in which X, Y, and Z are the axes of a right- 
handed, earth-centered, earth-fixed, rectangular Cartesian coordinate system 


z 



Figure A-1.- Geometrically determined longitude of a point on earth* s 
visible horizon relative to location of a satellite. 
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where the Z-axis is parallel with the earth's spin axis, the X-axis is parallel 
with the intersection of the Greenwich meridian and the equatorial planes, and 
the Y-axis completes a ri^t-handed system; Aq and Ay n are the longitudes 

of Ti and Vn, respectively, and are conventionally measured from the 
Greenwich meridian positive in the eastward direction; <f>^ is the geodetic lati- 
tude of Ti conventionally positive north from the equatorial plane; r^ and 
r^ are radial distances from the earth's spin axis to Tj^ and Vn^ respec- 
tively; the plane TiV^I is tangent to the ellipsoid at Tj^; and Nj is the 
normal to the ellipsoid at Tj^. 


From figure A-1 



(Al) 


^n,i \r,n “ ^i 


and 


"^i,n - ^,n " ^i 
+ (% sin 0i - Zi + 

From the law of cosines^ 


(A2) 

(A3) 


COS 


V T- 


- AZ 


1 




- 




m 


Thus, after substituting equations (Al) and (A2) into (A3) and the resulting 
expression into (A4), the result reduces to 


- Az 


cos AA- 


[Sin 0^' 
i^’^^cos 0^ 


n,i 


v,n 


(A5) 


But 


and, from figure A-1 


tan 0| = (l - e^)tan 0^ 





30 


APPENDIX A 


Therefore , 

sin 0j_ 

cos 0i ri(l - e2) 

and, after substituting into equation (A5) and factoring, 


cos AXn^i 


‘•v,n 


1 - 


^1-^1 ,n 
(l - e2)rfj 


From figure A-1, 


2 2 2 
rf=xf+ yf 


But, from the equation of an oblate spheroid, 


rf = X? + y? = ^(b^ - z?) 


(A6) 


(A?) 


Therefore, when equations (A 7 ) and (A2) are substituted in (a 6) and the results 
reduced 


cos = 


,2 

a ^ " ^v,n^i 

^ (b2 - 


(a8) 


It is convenient here to define two additional coordinate systems, 

X',Y',Z' and X" ,Y" ,Z" , which are, also, right-handed, rectangular, and 
cartesian. As shown in figure A-2, the X',Y',Z' system is earth-centered and 
the Z*-axis is parallel to the earth's spin axis and positive north, which is 
equivalent to the Z-axis of the X,Y,Z system. The Y'-axis lies in the equa- 
torial plane and parallel with the intersection of the vehicle's meridian and 
the equatorial plane and is positive in the direction of the satellite. The 
X'-axis completes the right-handed system. The X" ,Y" ,Z" coordinate system 
is vehicle -centered an^ defined in such a manner that the Y"-axis i_s parallel 
to the radius vector Rv,n positive in the same direction as Rv,n> 

Z"-axis lies in the plane of observation normal to the Y"-axis, positive in the 
positive direction of Dn^i and Pj_ degrees clockwise from the meridian of 
the vehicle; and the X"-axis completes the right-hand system. 

The coordinates of the satellite Xy y^ and Zy in the X,Y,Z 
system become 0, and Zy^n X',Y^,Z' system where 

^v,n “ ^,n ^v,n 


^v,n ~ ^,n ^v,n 
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Figure A-2,- An optical observation relative to a satellite- 
centered coordinate system. 


and the coordinates of the point of tangency Tj_ become y\_$ z-^. 

Further, since the vector D^^i (^n,i ~ ^n'^i originates at the 

satellite and terminates at the point Tj_ on the visible horizon, the coordi 
nates of Tj_ can be expressed in terms of D^^j_ and On,i 
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or 


y± = -Dn,i cos 

sin 

Thus, from figure A-2 

Pn,i = “^fi ■•■ (^v,n " ^i)^ '*' (^v,n “ 

®n,i “ ^n,i ’ ^n,i “ (^i) (^v,n " ^i) ■*" (^v,n 

= + Ny^n “ ^ (^v^n^i ’*' ^v,n^i) 


-Hf 


and, since 


. «2 1 / 2 
yl = r. cos AA^ = ^ — 

'^1 1 ^ v2 r. 


■b"^ ■"v,n 


(t>^ - Zv,n^i) 


equation (A9) "becomes 

^2 


^n,i = 4 + '^,n- 2a^ " ^Zy^n^i 1 - ^ 


(A9) 


(AlO) 


(All) 


Furthermore, by the cosine law, 

sS a. 1^2 


^,n \,i ■ 2Rv,n^,i cos - % 


or 


2^,n^n,i cos - IV,n ^n,i “ % 

O 

which, upon substituting equation (All) for i^,i and solving for 
reduces to 


^n,i = 


^v,n - - ^2^ 

Ry,n cos o^^j_ 


(A12) 


To find x^, yJ, and z^, first rotate the X" ,Y" ,Z" system about the 
Y"-axis through the angle (3j_, then rotate this system about the X"-axis through 
the angle 0^, and, finally, translate the origin to the earth* s center. In 


35 


APPENDIX A 


matrix notation 


r 


’ 




— 


— 


r 

~\ 


r 

4 


1 

0 

0 


cos 

0 

-sin 


0 



0 

4 

[ = 

0 

cos 0^^n 

-sin 0;^n 


0 

1 

0 


-Dn,i cos 

On,i 

■> * { 

^v,n 

zi 


0 

sin 0y,n 

cos 


sin 

0 

cos Pj^ 

1 

^n,i ®ln 



^v,n 

L > 


— 


— 




J 





V J 


and 


1 

Xi = 

“^n,i 

sin 

chi,i 

sin 

Pi 






! 

7i = 

-l^n,i' 

^cos 

“n,i 

cos 

0;,ni- 

sin 

Chl,i 

cos 

sin 0;^n) + ^v,n 

S (A13) 

^i = 

-Cn,i| 

cos 

°h,l 

sin 

0;,n - 

sin 

°^,i 

cos Pj^ 

0v,n) + 2y,n_^ 



Substituting equation (A12) for ^ in the expression for and 

solving for Zj_ yields 

- — - III ! -1 (A14) 

(sin^0v^n ' ^,i Pi K ,n ^os (j>'^ 

which is the desired equation. 

From figure A-2^ 

X % = pRv,n% ^l,n 

where p is a unit vector normal to the plane formed by Ry^n ^i" Thus, 

(^,n X %) • (Rv,n X %) 

1 = p • p = 

^^,n^i ^i,n^ 
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which, after substituting the rectangular representations of 
perfonning the indicated operations, and solving for sin^rij^ becomes 


. 2 ^ 
sin = — 


,n4 - 2^v,nS,ny^i + 4 ,n + ^ 4,nH)‘ 


Now, by the sine law. 


44, n 


2 2 

o R,- sin Ti . ^ 

sin^c. . = ^ 

n,i 


Thus, 


. 2 V 

sin = 


,A - * {y'±f4,r, * ^v,n(^i)^ + 4,nAy 




Substituting the identities 


^v,n = IV,n sin 
^v,n = IV,n cos 


I _ a^ 1 A 2 \ 

^i = ^ - ^v,n^i) 


* i /2 { « \2 

Xi - \^y^) = 


a"^/, 2 2 \ / f \2 

- HJ - (^i) 


1/2 


(1 - f)^ 


1 _ ^ 

b 2 (1 - f )2 


and rearranging the terms gives 
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is obtained. The desired expression for is, therefore. 


^-l[ ~^n,i /^,i ~ ^^n,i‘^n,i I 

\ ^n,i j 


°4i,i = 

The positive sign is chosen so that On^i is a minimum when Pn^i is 


(A17) 

zero. 
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AN ERROR EQUATION FOR a 


In the section entitled "Assumptions and General Theory/' the general 
expression for the error equation of alpha (eq. (l 5 )) was derived, and in 
appendix A, the expression for F was derived. To coiiiplete the solution of 
equation (l 5 )^ it is necessary to derive first the expressions for the func- 
tions A, and P so that the proper partial derivatives for determining 

the coefficients of equation (ll) can be obtained. 

The partial derivatives of 0, A, and P are unwieldy when no form of 
generalization is utilized. In an effort to generalize, the necessary partial 
differentiation of individual expressions will be indicated only in this sec- 
tion. The reader is, therefore, referred to appendix C for the necessary alge- 
braic expressions which result from the partial differentiation. 

Before A, and P can be derived, it is necessary to decide what 

method will be used to fix the position of the satellite. Because of the high 
accuracy obtainable, radar was chosen for this paper. 

There are several techniques for fixing the instantaneous position of a 
space vehicle by radar observations in use today. Since it is impractical to 
discuss each technique and its individual application to the theoiy of this 
work, the conventional measurements of range Rsv,n^ azimuth Azg^j^, and the 

elevation above the horizon 7 ^^ are taken to be observables. 

Of the three observables listed, only the range is independent of the 
deflection of the vertical. The dependence of ^n compo- 

nents of the deflection of the vertical will be shown subsequently in this 
section. 

The relations of equation ( 6 ) can be rewritten as follows: 




+ ^ dr„ t — dAz. 




o7n ^^^sv,n ^ 


dA- 


da + — df + ^ 

oa hf S0s 


•v,n = — 'ia + ^ df + ^ d0s + ^ dAs + dhs + <i7n + ^ 








5hg u/n ^^sv,n 


(Bl) 


■ I - * I “ * I; ^ 1 * It * I; ^ 


5A.q 


3hp 


^7n 
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The partial derivatives of A, and P with respect to the range are zero, 

as explained previously in the text. 


Thus, the rectangular coordinates of the radar station in the X,Y,Z 
coordinate system are 


and 


Xg = (Ng + hg)cos 0s sin Xg 
yg = (Ng + hg)cos 0g cos Ag 






Ns 


Ws 


(B2) 


(B5) 


is the expression for the magnitude of the normal to the ellipsoid vhere 


w| = 1 - (pf - f^)sin^0£ 


(ref. l) . Therefore, 
tan 0g = 


^x2 + y2 


1 - 


N 


® (2f - f^] 


(Ns + hg) 


tan 0£ 


(b4) 


where 


2f - f^ = 


is the relation between the geocentric and geodetic latitudes. 
The total differential of equation (B3) is 




SNg 

^s 


df + d0s 


and that of (b 4) is 

^ <iNs ah, . ^ df * 


ajzli 30; 30; 30; 


She 


Sf 


d0j 


(B5) 


(B6) 


Thus, after substituting equation (B5) into (b 6) and collecting coefficients of 
like terms 
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SNs ^ 


da + 



(B7) 


In figure B-1, the^ relation at the radar station (or site) s of the nor- 
mal to the ellipsoid. Ng, the geocentric radius vector of the site Rs^ and the 
range vector Rsv^n is shown. (The arrowheads indicate positive directions of 
the vectors.) The vectors Ng and Rgv^n plane of the paper so 

^hat the line segment sH represents the plane of the horizon and the vectors 
Ng and Rg are in the meridian plane. 


If s is considered to he the center of a unit sphere, the spherical 
triangle DEF in figure B-1 can he formed where the angle at D is 


D = n - Azg^n 

the angle at E is the spherical azimuth of the satellite Azg^n^ 
given hy the relation 

S0s = 0s - 0s 

Thus, from the cosine law for spherical triangles, 

sin 7^ = cos 60g sin 7^ - sin 50g cos 7^ cos A.Zg^^ (B9) 


and from the law of sines 


cos 7n 

sin Aza „ = sin Azg^n 


cos 7. 


n 


The total differential of equation (B9) is 

o7n ^Zg^n °(S0g) 

and that of (BIO) is 




dAz 


SAz. 


dAz 


s ,n 


3iAz 


2“ dAz + dr' + d-v 

^2g,n + °-'n 


g^n 


n 


^7n 


(BIO) 


(Bll) 


(B12) 
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As stated previously, Azg^n independent of the deflec 

tion of the vertical. In reference 1 it is shown that 


dA.Zg^n “ Az^ 


(B] 


The egiiivalent expression for must he derived. 
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The normal to the best fitting ellipsoid will be deflected At] sec of arc 
in the meridian of the normal to the old ellipsoid and sec of arc in the 

prime vertical. It is now desired to determine the change in the elevation 
d7j^ which will result. 

Since the variables involved are arcs and will be measured at the common 
origin s (the radar station)^ figure B-2 represents the geometry of the prob- 
lem to be solved. The angle included by the arcs At) and A| is a ri^t 
angle. Thus, since At] and A| are small, the change in the deflection of 
the vertical v can be written as 

= A^2 + At)^ 



Figure B-2.- Effect on elevation resulting from a change in deflec- 
tion of the vertical. 
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Further, the law of cosines gives 


cos 


I - (7n + ^^n)l = " :>'n)cos v + sin v sin^| - 7n)cos(Azg^n “ ^^) 


which reduces to 


sin^7n ‘^^n) ~ ^n ^ ^ 7n(*^°® -^^g,n M- + sin Azg^n 1^) 


where 


li = tan"^ ^ = sin"^ ^ = cos"^ ^ 
At) V V 


In reference 1, A^ and At] are related to the corresponding changes in 
latitude and longitude by the expressions 

A| = dXs cos 0s 

At^ = d0s 

(in ref. A| and At] are defined to be positive in directions opposite to 
those used here. Thus^ the signs in the above expressions must differ from 
those in the reference.) Thus, since v is small and dy^ will be small, 

d7n = <^os ^^s ^^g,n ^s (B15) 


For expediency, the coefficients of d0g and dXs in eq^uations (B14) 
and (BI 5 ) are defined as 


clAz 


sin Azg^n tan 7n = ^ 




(BI 6 ) 


cos Az 


_ ^ 
= S03 


SAz, 


g ^ 

sin 0s - cos Azg^n. 0s 7n = ’ 


Thus, 


/ ^7n 

sin Azg^n cos Ps - 


SA2 


dAz 


'g,n 


g^n 


50, 


d0, 


5Az 


g^n 


5Ac 


dAc 


(BIT) 
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and 




By definition 


50s = 0s - 0s 


Thus , 


d(S0s) = <i0s + i0; = d0s - i0s 


and, after substituting equation (B7) for 


<i(80s) = -^da-^df-^dhs + |l 




She 


^sV 


d0£ 


(I 


Finally, substituting equations (BI 9 ), (BI 8 ), and (BI 7 ) into (Bll) and 
(BI 3 ) £ind collecting coefficients of like terms give 


t ^^n ^^n ^^n a ^^n ^^n 


Sf 


S0C 


SAc 


She 


where 
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where 


^Zg,n _ 

^2g,n 

^k 

% 

Sa 

^yk 

S(S0s) 

Sa 

^^s,n _ 

^^s,n 

S 7 ' 


Sf 

K 

S(6j^g) 

Sf 

SAzg,n __ 


K 


Shg 

K 

S(S0s) 

Shs 

_ 1 



n^' 

^S ' 

V ^rn 




^ ^^s,n 

K 

L ■ 




If-; 

II 

•N 

CO 


^^S,l 


SXg 


^yk 









bkz. 




\bkz 


(B23) 


g^n 




ihkz 


g^n 


SAc 




The rectangular coordinates of the radar station are given in terms of the 
geodetic latitude, longitude, and station height by equation (B2) . The rec- 
tangular coordinates are given as functions of spherical coordinates by 


Xs 

“ ^s 

cos 0s 

sin 

> 

As 


Xs 

= Rs 

cos 0s 

cos 

^s 

) (B24) 

Zs 

= Rs 

sin 0s 


J 


1 of 

equations 

(B2) 

and (B2^) are equated, the 


resulting expression solved for R 3 would be 

/ V cos 0 


(B25) 


The same expression would result with the coordinate. The Zg coordinate 

would, however, give a different expression for Rg, but the numerical value 
would be the same . 

The total differential of Rg is 


Sr 


Sr 


Sr, 


Sr, 


^ dUg t ^ dhg t ^ d0g ^ ^ 0 ; 


Sh 


S 0 £ 


s0; 


(B26) 
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Since dNs and d0g are given by equations (B5) and (b 6), respectively, after 
substituting into equation (B26) and collecting the coefficients of like terms, 


/ap dR a0;\dN 

dRc; = — ^ ^ ^ da + 

IdNs d0: dNs/da 


S0s 301 30s 


d0£ 





dhs + 


W dNgi^ 


dNc 


^^s ^s\aNc 


dR d0' 
s 

d0; dNg/df ■ d0; df 


df 


(B27) 


Consider figure B-3 in which the polar triangle (NP)(Vn^g )(s) is formed 
by the intersection of the meridian planes of the radar station and vehicle and 
the plane containing the geocentric radii of the vehicle and the radar station 
with a sphere of radius Rg. The angle at NP is 



Figure B-5*- Spherical coordinates of a satellite relative to a 

radar station . 
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^sv,n “ \-,n ~ 

and that at s is the spherical azimuth Azg^n satellite. Further^ 

the arcs NPs and NPV^ are 



and 


NPV. 


= 2L - 

n 2 *^v,n 


The arc sV^ is defined as 9n and^ from the plane triangle of fig- 
ure B-5, it is seen that 


tan 0n = 


^^sv,n 7n 
Bs Bsv^n ^n 


(B28) 


Since the radar observation of range is considered herein to be errorless , the 
total differential of 0n is 


80 80 


(B29) 


But, d7^ is given by equation (B20) and dRs is given by (B27). Thus, after 
substituting and collecting terms 


80 „ 80„ 80„ ^ 80^ 50 „ 


(B30) 


Wl? 0"*^^ 


^ ^ ^ 80;\8n3 87A B0^ 

8a 8Rsy8Ng 8Nsy8a 8a 8^ 


n 


8f 

^Qn ^ 
^s 



_ 

i' 


8Ra 8f 87 ^ / 


^Ba . ^B^ 801 


8 Rs 80 s ^ 7 n 


(B31) 


(Equation continued on next page) 
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SAs SAs ^7n 

^ ^ + 5zi 

SNs ^Rs\5hs ^s ^hg/ Stis ^7n 




(B51) 


Tvfo additional expressions which can be obtained from the spherical tri- 
angle of figure B-5 are 


and 


sin 0v^n ®n 0s ®n cos 0g cos Azg^n 


sin Az 

sin AAsv,n “ ' ■ ^ sin 9n 


cos 0 ^ 

'^Y,n 

The total differential of equation (B 52 ) is 

flfA' - -t- ^v,n ,Q ^0y,n 

^K,n - <i0s + <ien + _ *^^s,n 

and that of (B55) is 


(B52) 


(B35) 


aen 


(B54) 


s ,n 


dAA^ 


sv.n 


ClAAg„ „ n J| n 

d0n + - V dC . + . ^ . . Jy .2R dAz. 


Se 


n 




<,n 


ciAz 


(B55) 


s,n 


The differentials d0g, d9n^ and dAzg^j^ are given, respectively, by 

equations (B7), (BJO), and (B22). Therefore, after substituting into equa- 
tion (B3^) and collecting coefficients of like terms. 


<i0v,n = '^a,n ^ + ^0g,n ^0s + <^Ag,n + <^hs,n 


(B36) 


where 




t 

= 



8n„ 

— t 

30;, n 

^n ^ 

^v,n 

^2g,n 


d> ^ XI. 


8Ns 

Sa 

89n 

Sa 

aAzg,n 

bsi 


,» 

^v,n j 

f^'s 

^Ng 

^s\ 

^ 30;,n ae„ ^ 

30;,n 

aAZg 

f,n 

■ 

^SNg 


3f / 

d9n 

8f 

^^s,n 

8f 


. KJ^'s 

^ J 


.^4 

,n ^n 


8Az, 

0s>^^ 


\8Ns 

. a0g 


^n ^n 

^^s,n ^0 




(B37) 


(Equation continued on next page) 
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_ ^v,n ^ ^v,n 
SAs ^Azg^n 

^0 ’ S0 * S0 * ^0 ^0 * hkz 

^s _j_ ^Y,n Q^n ^ ^v^n s,n 

S0g Shg ^©n SA.Zs^n 




(B37) 


Before the functional coefficients of the Ay^n error equation can be 
derived, it is necessary to refer to the definition of ^!i^sv,n" By definition 


Thus, 


and 


^sv,n “ \r,n " ^s 


^^^sv,n “ '^^v,n “ <^^s 


(B38) 


dAy,n - '^^sv,n *^s 


(B39) 


After substituting equation (B35) into (B39)> By a procedure similar to 
that for deriving the coefficients of d0y^j^. 


dA. 


v,n = K,n ^ + ^,n (Bi^•0) 


where 


Aa,n ~ ^a,n 




sv,n 


SiAA 


sv.n '^''n 




hkz^ 




sv.n ^ "s,n 




80 n 8a 8 Azc 


8a 


A» _ a' ^^sv,n ^*^sv,n ^®n ^^sv,n ^^s,n 
^,n - f,n SAzs^n 

Aa - ^^sv,n ^^sv,n ^^sv,n ^^s,n 

0s ^v,n 8Azg^n 80 g 


As^n ^ Ag,n 8Azg^n 8Ag 


8^^ 


+ 


8^ 


86_ 8aA, 


sv_^ n ^ 


8Az, 


sv,n s,n 


Ai 


Bo^n 


= 


8^ 


sv,n 




80V, 8;:^ 


8Az, 


sv,n '-^n ^ sv,n s ,n 


80|^j^ 80n 8hg 8Azg^n 8hg 


(B4l) 


To complete the derivation of the coefficients of equation (8), an expres- 
sion is needed for dRy^n* From the plane triangle of figure B-3 
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^ ^ cos 7 ^ 

- Rsv,n gin 0n 


sin 9n = 


_ cos sin 


sin Azc 


Thus, 


^,n ~ ^sv,n 


cos 7' sin Az^ 


K,n 


sin AXc 


The differential of equation (b 43) is 


where 


p' - ,T,' 4. A' + ^\,n ^\,n ^^s,n 

a,n a,n^,^^ ^ ^ 7 ^ Sa ^z^^n 

^sv,n ^^s,n 

pi = $1 + A» ^7n ^ 

0s S0,V n ^s^’^ Sk^sXgy^n ^7^ ^s ^ 


's,n ^s 


pi „ = $1 + /"a* _ i'\_^LdL + ^ ^^s,n 

^A,n ^ /^sv,n S7A ^Zg^n ^^s 


^v,n ^^^sv,n ^7jj Shg ^Azg^n ^^s 


(B43) 


^ * fe: ^ <,n - ^3v,n (B^>*) 


After substituting for ‘^Zg^n, ‘^0v,n^ *^7n^ 

■®v,n = K,n ^ * Bf,n + B0„n ■*»*= * ^{^,n '^s + P;^,„ (B“5) 


Now that the coefficients for equations (8) have been determined, it 
remains only to determine the coefficients of (12). It will be recalled that 
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where ®n,l^ ^n,i given hy equation (A16). As a result, 

^ ^,i ^V,i ^ ^,i ^^n,i ^ ^,l ^n,i 

\^/n,i ^n,i ^ ^®n,i ^ ^n,i ^ 

^ ^,i ^n,i ^ ^,i ^ ^,i ^n,i 

^f/n^i ^n,i ^®n,i ^n,i 

bF \ ^ ^,i ^,i 

^v,n/n,i ^^n,i ^v,n 

/ ap \ ^ ^,i ^,i 

\^^^v,n/jj^j^ ^n,i ^^,n 

which are the desired expressions. 

The error equation for oc^ j_ is given in a general form by equation (l3)* 

As stated previously, the coefficients of equation (15) are functions of those 
of (8), (9), and (12). 

In this section the coefficients of equations (8) and (l2) were derived 
as functions of indicated partial derivatives and are given in equations (B37 )^ 
(b4i), (b 46), and (b 47). The indicated operations have been executed and the 
results are given in appendix C. To obtain the algebraic expressions for 
equations (B37 )) (b 46), and (B^7), it is required only to make neces- 

sary substitutions. 


^,i 

aBn,l ^ 

8°4i,i 

8Cn 

A 

8®n,i 



90V 

A 

8on,i 




A 

8Bn,i 

^^v,n 

8*^n,i 

8Rv 






J 
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ANALYTICAL EXPRESSIONS FOR VARIOUS PARTIAL DERIVATIVES 


In the section entitled "The Error Equation for a" the coefficients of 
equations (8) and ( 15 ) were derived as functions of partial derivatives. In 
this appendix the partial derivatives of those expressions on which the coef- 
ficients are dependent will be given. 


The partial derivatives of expressions (B3), (B9)? (BIO), (B25)^ 

(B28), (B 32 ), (B33)^ -'ind (b 43) are as follows: 


80S 



s 


8Ns _ Ng 

8a a 


8Ng 


= (1 


f)^^| sin20s 
W2 


8hg 

8Ns 


= (2f - f2) 

(Ns + 

= -(2f - f2) i!; 

(Ns + 



tan 0g cos20g 
• tan 0g cos^0 


I 

s 


80s 


1 


N, 


No + h. 



I 

s 


s 


80^ 


8f 


- 2(1 - f) 


Ns 

Ns "*■ ^s 


tan 0£ 


cos^0 


I 

s 


8Rs ^ Rg 

8hg Ng + hg 

8Rs ^ Rg 
8Ng Ng + hg 



Rg tan 0g 


52 


APPENDIX C 


SRs 


-Rg tan 03 


Se 


n 


-sin^Gn 


SRg Rsv,n 


^0n 

^7* 

>-’/n 


Ra \ sin^9„ 


Sin 7n + ^]—2l 
\5'sv,n /cos^7n 


sin 60g sin 7^ + cos 60g cos 7^ cos Azg^n 
S50S cos 7n 

cos 60g cos 7n + sin 60s sin 7^ cos AZg ^ 

^7n cos 7^ 

sin 60s cos 7n sin Azg^n 


SiAz 


g^n 


cos 7 n 


SAz 


s»n ^ , . » 

= tan Azg^n ^n 




6 Az 


^7n 

hkz 


= -tan Az„ „ tan 7. 


s,n 'n 


cAj 


— = AZg^n AZg^n 

g^n 


^S^v,n ®n ^s “ ®n 0s -^^s,n 

~b^ cos 0^”n 


0s ®n + cos 0g cos tn cos Azg^^ 
__ cos 0;^j^ 

^v,n ^ -sin Gn cos 0^ sin AZg^^ 

^Zs^n cos 0y^n 


55 


APPENDIX C 




sv^n 


Se 


- tan cot 0^^ 


n 




SAz 


= tan AAg^ cot Az 

s ,n 


s ,n 




^sv,n "tan 

^Y,n 


_ _ -tan >' 

'n 

_ . R,,, cot Az,,, 
_ 


S0. 


7^ = Rv,n tan 0 ;,, 


v,n 


,n 

IV,n =°‘ ^sv,n 

'-*-^sv,n 

The partial derivatives given in equations (BI7) and (BI8) are 


= sin Azg,n tan 7 n 


dAz 


S^n 


= sin 0s - cos Azg,, cos 0g tan 7, 


^^n A 

^3 ' 


= Sin Azg^n cos 0g 
'S 


If the partial derivatives of ciji i respect to ^^i^ 

^n,i obtained and substituted into expressions (b 47 )> the results will 

become 
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— = - ( 2 f - f 2 )sinVv,n ' & ' 2 ( 2 f - - -|i-yin 20 ;^^ 

\ V ^V,n/ \ ^v,n/ 

+ |& - 2(2f - f^O^l - ^j^cos20;^^ cos^p^^ 


- <(2f - f2)(l - 


Sa, 


n,l ^ 

SRv,n 




^.nj 


Vl - f)‘^ + (2f - f2)[(2f - f2)sin^0^^n + 2(1 - 


+ (2f - f2) (2f - t^)sin%^^ + {1 - fffl - cos20;^jj cos2p^ 

L ’ \ ^,n/J I 




^ = 2R?,n(2f - f2)/(2f - f2)sin20;^n + (l - f)^(l - + | 

'v,n V Rv,n/ 


^,n 


(2f - f2)cos20;^n 


- (1 - f) 1 


^,n. 


cos^pAsin 


^ = 2a(2f - 




SBn,i p 
__^ = 2 a‘=^' 
df 


f2) (2f - f2)sin20;^n + /l - 2 -|^Vl - f)‘ 


sin 20y^n cos p^^ 


(1 - f)<|2(2f - f2)sin20;^n + [l - 2(2f - f2)] ^1 - ^Lj^sin 20^^„ cos pi 


= 2(1 - tf {2f - f2)-^ sin 20; cos p^ 
^^y,n r3 




_Ih± = a2^(2f - f2)sin^20;^n + 2 

^v,n 


(2f - f2)sln20;^^ + (1 - f)^^l - -|^j cos 20;^^cos Pj^(2f - f^) 
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= -2a// 1 - 2 


Sa 


R|,n, 


(1 - ff + 


2 1 - 


r 2 


(1 - f)^ + (2f - f2)sin20^^^ 


(2f - f2)sin^0^ 


^n,i 

~W~ 


2a^(l - f) 2(1 



f)^ 




+ 2(2f - f2)sin20;^Asin20;^j^ 


^n,i 


.2 ^^{1 - ff\{l - + (2f - f^)sin^0^^n 








= -a2(2f - f^) 


(1 - f) 2 - 


r 2 


2(2f - f^)sin^ji 


v,n 


sin 20^ 


v,n 


Finally^ it can be stated that since 

^n,i = ^0v,n^ ^,n) 

and expressions for d0 and dE are given by equation (lO) as fiinctions of 
those relations given in (B57)# (b 46), the coefficients of (l5) can 

be readily obtained. 
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COEFFICIENTS OF THE DE GRAAFF -HUNTER EQUATIONS 


The method of de Graaff -Hunter was chosen over that of Vening Meinesz 
since it is simpler to work with. This method is given in reference 1 on 
pages I 27 -I 3 O in the form most widely used. 

The accuracy of the two methods is largely dependent upon distance from 
the origin of the spheroid. For shorter distances the two methods will give 
accuracies of 0.001" of arc. For greater distances, the accuracy of the method 
of de Graaff -Hunter falls off whereas that of Vening Meinesz will give accura- 
cies of 0.001" of arc or better. 


The form of the de Graaff -Hunter equations as given in reference 1 is not 
satisfactory for this paper. The form of equation ( 9 ) is that form desired. 

To obtain it, the equations were expanded, and the corresponding coefficients 
of the five arbitrarily defined constants of a datum were collected. In so 
doing, it was determined that by making certain trigonometric substitutions the 
resulting coefficients can be simplified. 

Further, in reference 1 SNg and 6 Nq were considered positive when the 
new spheroid was above the old. In this paper hg is positive for elevations 
above the spheroid. Since hg and 5Ng increase in opposite directions, 

Ahg = -A5Ng 


Thus, the expressions of the de Graaff -Hunter equations giving the change in 
spheroidal height must be multiplied by a minus one. 

The coefficients of equation ( 9 ) are as follows: 

_ (cos Ug sin Uq - sin Ug cos Uq cos cd)c 


Lf = c 


sin 2ug - sin Ug cos Uq sin^Uo cos cu - cos Ug sin Uq^I + cos^uiq^ 


= c sin Ug sin co 


L| = -c^sin Ug sin Uq cos co + cos Ug cos Uq^ 


"6N 


c(sin Ug cos Uq cos cd - cos Ug sin Uq) 
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Ma. = 


Mf = 


Mr, = 


= 


yg COS Uq sin 00 


-^1 - sec 0s cos Uq sin^Q sin co 
-^1 - ^^sec 0s cos (D 

■(^ " si’ 

(l - ^)sec 03 


sin Uq sin ou 


cos Uq sin O) 


Sg, = cos Ug cos 0 ) cos Uq + sin Ug sin Uq - 1 

Sf = -ajsin Ug sin Uq + cos^q) - cos Ug cos a> sin^UQ cos Uq - sin^ 


= -a cos Ug sin co 


S| = -a(sin Ug cos Ug - cos Ug sin Uq cos co) 


Sgjj = -(cos Ug cos (n cos Uq + sin Ug sin cDq) 


where 

c = ^1 - ^) (l + f cos20g) 
tan Ug = (l - f)tan 0g 
tan Uq = (l - f)tan 0g 

and 


a> - As - Aq 
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